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How common is large-scale entanglement in nature? As 
a first step towards addressing this question, we study the 
robustness of multi-party entanglement under local decoher- 
ence, modeled by partially depolarizing channels acting in- 
dependently on each subsystem. Surprisingly, we find that 
n-qubit GHZ entanglement can stand more than 55 % local 
depolarization in the limit n — > oo, and that GHZ states are 
more robust than other generic states of 3 and 4 qubits. We 
also study spin-squeezed states in the limit n — > oo and find 
that they too can stand considerable local depolarization. 

Entanglement is certainly one of the most dramati- 
cally non-classical features of quantum physics. There 
are states of composite quantum systems that cannot be 
decomposed into probabilistic combinations of product 
states: there is no way of writing the density matrix of 
such a state p in the form p = J2Pi a l ® °2 ® ■■■ ® a n< 

i 

where a], is a state of the fc-th subsystem and the Pi > 
are probabilities. Such states p are called inseparable or 
entangled. One could say that for entangled systems the 
whole is indeed more than the sum of its parts. Only en- 
tangled states can exhibit quantum non-locality pL In 
this case there is no way of reproducing the predictions 
of quantum physics with classical systems, unless there is 
instantaneous communication over arbitrary distances. 

Entanglement has recently been studied extensively in 
the context of quantum computation and quantum com- 
munication Q . Quantum computation would require the 
creation and maintenance of highly complex entangled 
states of many subsystems. This is a difficult task be- 
cause of decoherence [9 . A system interacts with its sur- 
roundings, which creates entanglement between system 
and environment, at the same time reducing the entan- 
glement within the system itself. 

Such considerations lead to a natural physical ques- 
tion: how common is complex entanglement - entangle- 
ment between many subsystems - in nature? Are there 
natural systems that contain substantial large-scale en- 
tanglement between their constituents? There seem to 
be two conflicting intuitions in response to this question. 
On the one hand, entanglement seems to be very fragile 
under decoherence, which is emphasized by the difficulty 
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of building a quantum computer. On the other hand, 
measure-theoretic considerations suggest that the set of 
inseparable states is much larger than the set of separable 
ones |Q , which seems to imply that entanglement should 
be relatively common. 

One might hope to gain some insight into the above 
question by analyzing the behaviour of complex entan- 
gled states under decoherence. in particular by studying 
how fast the entanglement disappears. This is not an 
easy task, because it requires practical criteria for insep- 
arability. However, for some special sets of states such 
criteria have been found (^||, which we will use in the 
following. 

We will consider the evolution of certain multi-qubit 
states under the action of partially depolarizing chan- 
nels, to be defined below, acting independently on every 
qubit. The multi-qubit states are our models for com- 
plex entanglement, and the depolarizing channel is our 
toy model for decoherence, corresponding to local and 
homogeneous noise. We ask how much local depolariza- 
tion is possible such that the states are still entangled, in 
this way quantifying their robustness under decoherence. 

Our notion of robustness is thus different from the one 
of , who considered global mixing of entangled states 
with locally prepared noise. Admixing locally prepared 
noise is different from local decoherence; it corresponds to 
a scenario where the local "decohering agents" commu- 
nicate when to add noise, and when to leave the system 
undisturbed. 

It also differs from the notion of entanglement persis- 
tency proposed in J8|, which is defined as the minimum 
number of local measurements needed in order to defi- 
nitely disentangle a given state. This definition involves 
an optimization over all possible local measurements. It 
thus corresponds to intelligent and informed adversaries 
who optimize coordinated decohering measurements in 
order to most efficiently destroy the entanglement of a 
state. The depolarizing channels which we consider here, 
applied to each subsystem, correspond to independent 
measurements in random bases on each local component 
of the entangled state. 

We will consider qubits with basis states |0) and |1). 
The completely depolarizing channel is defined by the 
transformation — > <%|l, where 1 = |0)(0| + |1)(1|. 
When applied to a subsystem A of a composite sys- 
tem AB this corresponds to the transformation pab — > 
\ \ ® Tr a Pab- All correlations between the system and 
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its environment are destroyed, and the system itself is 
put into the completely mixed state. 

The partially depolarizing channel Cd is defined by ap- 
plying the completely depolarizing channel with a proba- 
bility d, and applying the identity transformation with a 
probability 1 — d. This corresponds to the transformation 
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If the density matrix of the system is written in the basis 
of Pauli matrices, one can easily convince oneself that 
(|l|) corresponds to multiplying all <7j by the scaling factor 
s = 1 — d, while the unit matrix 1 remains unaffected. 

An interesting physical realization of a partially depo- 
larizing channel is by random measurements. One can 
show that if a qubit is subjected to a Von Neumann 
measurement in a basis chosen uniformly at random, the 
effect on the density matrix corresponds to a partially 
depolarizing channel with d = 2/3. 

We will refer to d as the amount of local depolariza- 
tion. Denote the product of partial depolarizing channels 
with depolarization d on each qubit by Cf n . For a given 
n-party entangled state p, we will be interested in the 
critical amount of depolarization d cr n where Cf n (p) be- 
comes separable, or - in the absence of a necessary and 
sufficient condition - ceases to fulfill certain sufficient con- 
ditions for entanglement [||J^] . In both cases d cr a (ft) will 
be a quantitative signature of the robustness of the en- 
tanglement in that state. 

There are two families of states where we have been 
able to obtain explicit results in the limit of n, the num- 
ber of subsystems, going to infinity, namely GHZ states 
and spin-squeezed states JlO| of n qubits. Our main re- 
sult is that, surprisingly, in both cases a very substan- 
tial amount of local depolarization does not destroy the 
multiparty entanglement. We have also compared GHZ 
states to other generic entangled states for n = 3 and 
n = 4 qubits, with (to us) counter-intuitive results. 

For the following discussion we find it convenient to 
define P = |0> <0] , Pi = |1)(1|, a+ = |0)(1|, and cr_ = 
|1)(0|. In this basis the partially depolarizing channel (|l|) 
corresponds to the following transformations: 
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Recall that s = 1 — d. 

Let us first discuss the case of GHZ states. An n-qubit 

GHZ state l/\/2(|00 . . . 0) + |11 1>) has the density 

matrix 



p = I(P® " + Pf n + <jf n + (7?"). 
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Application of the channel (g) to every qubit multiplies 
the off-diagonal terms af n and cr®" by s n . The diagonal 



terms P and Pf n give rise to new diagonal terms of 

and all permutations thereof, 



the form P, 
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for k ranging from to n. The coefficients of these 
terms are given by 
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In order to analyze the entanglement of states of the 
above form it is sufficient to study the entanglement for 
bipartite cuts || , where some qubits constitute one sub- 
system and all the other qubits the other subsystem. For 
each such cut one considers the partial transposition || . 
As long as the partial transpose of a state has negative 
eigenvalues, the state is definitely entangled. 

Take partial transposition of the first k qubits as an 
example. The diagonal terms do not change, while 
cr®" = (|0)(1|)®" goes into an off-diagonal term between 
the states |l)® fe |0)®( n - fc ) and |0)® fe |l)®("~ fc ). It is then 
easy to see that the state will have non-positive partial 
transpose for some bipartite cut, and thus definitely be 
entangled, as long as 
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for some k. In || it was moreover shown that if the state 
has positive partial transpose with respect to all bipartite 
cuts, then it is separable, so that the above condition for 
inseparability is both sufficient and necessary. 

The smallest eigenvalue A& is the one for k — n/2, if 
n is even, or k = m, for n = 2m + 1. Therefore we can 
immediately make a statement about entanglement after 
local depolarization for general n (consider even n for 
simplicity). The state is definitely entangled as long as 
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The critical value of s where the entanglement disappears 
is given by 
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for even n. It is very remarkable that s cr it decreases with 
n. The entanglement in the GHZ state thus becomes 
more robust under local depolarization when the number 
of subsystems is increased. In the limit n — * oo, s^ it = 
= 0.447. This means that a very large GHZ state can 
stand more than 55 % depolarization (d^ it = 1 — s^ it = 
0.553) of every qubit before it becomes separable. 

A natural question to ask is how the robustness of 
the GHZ state compares to other multiparty entangled 
states. We do not know of any other family of multi-qubit 
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n = 3 


1-2 cut 


G 3 


0.443 


w 3 


0.425 


n = 4 


1-3 cut 


2-2 cut 


G 4 


0.423 


0.489 


W4, 


0.423 


0.423 


X 4 


0.416 


0.453 


B 4 


0.468 


0.450 



TABLE I. Critical values of local depolarization d for 
which several generic entangled states of 3 and 4 qubits 
become PPT (with positive partial transpose). For the 
4-qubit states the values for both the 1-3 and the 2-2 bi- 
partite cuts are given. One sees that G4, the 4-qubit 
GHZ state, stays NPT (with negative partial transpose, and 
thus entangled) up to the largest value of depolarization, 
dcrit = 0.489. For the GHZ case PPT-ness and separabil- 
ity are equivalent. The states considered were the follow- 
ing: |G 3 ) = ^(1000) + |111»,|G 4 > = 75(10000) + |1111», 
\W 3 ) - ^(|001) + |010) + |100)), 

|w 4 ) = 5(10001) + 0010) + 10100) + 1 1000)), 

\X 4 ) = ^(|0011) + |0101) + |0110) + |1001) + |1010) + |1100)), 
and \Ba) = |(|0000) + |0011) + 1100) - |1111>). 



states for which there is a necessary and sufficient con- 
dition for entanglement. As a first step, we have there- 
fore chosen to apply the partial transposition criterion 
H to some characteristic states of 3 and 4 qubits. In the 
4-qubit case we have studied both the 1-3 and the 2-2 
bipartite cuts. These results can be directly compared 
to the results in the GHZ case, where the partial trans- 
position criterion happens to be necessary and sufficient. 
Our results are collected in Table |. The main - again 
surprising - conclusion is that GHZ states are compara- 
tively robust. 

The states considered were the following. For n = 3, 
the state |W 3 ) = 4-(|001) + |010) + |100)). For a symmet- 
ric state of 3 qubits, there is just one possible bipartite 
cut (1-2). The state |W 3 ) becomes PPT (with positive 
partial transpose) across this cut for a lower value of de- 
polarization than the 3-qubit GHZ state, i.e., at least 
according to the partial transposition criterion, it is less 
robust under local depolarization. 

For n — 4, we have studied three different states, 
namely \W 4 ) = |(|0001) + |0010) + |0100) + |1000>), 
|X 4 ) = ^g(|0011) + 10101) + |0110) + |1001) + |1010) + 

|1100)), and \B 4 ) = |(|0000) + |0011) + |1100) - |1111>) 
||. There is only a single instance where one of them 
is more robust than the 4-qubit GHZ state: the 1-3 cut 
for \B 4 ) remains NPT for higher values of depolarization. 
However, if the 2-2 cut is also taken into account, again 
the GHZ state remains NPT for higher d values. 

It is worth noting that the state \Xi), which can also 
be written as \X 4 ) = -i(|00)|ll) +21^+)^+) + |H)|00» 



originally has more entanglement across the 2-2 cut than 
the GHZ state: the Von Neumann entropy of the re- 
duced state is 1.252 compared to 1 for the GHZ state 
(for pure states, the Von Neumann entropy is a good 
measure of entanglement, cf. Q). Nevertheless it reaches 
the boundary of PPT states before the GHZ state. A 
more trivial example for the fact that the amount of en- 
tanglement and its robustness do not have a simple re- 
lation is provided by any number of independent shared 
2-qubit singlet states. Then the collective state is exactly 
as robust (or fragile) as an individual singlet, which has 
d C rit = l-w(2) = 0.423. 

Remarkably, there is also no positive correlation be- 
tween the Schmidt rank of the states with respect to the 
2-2 cut (i.e. the number of terms in the corresponding 
Schmidt decomposition) and their robustness under local 
depolarization, cf. |ll). The state \X 4 ) has Schmidt rank 
3, nevertheless it is less robust than the GHZ state, which 
has Schmidt rank 2. Again, an even simpler example is 
provided by a state of two independent shared singlets, 
\S 4 ) = |(|00)|00) + |01)|01) + |10)|10> + |11> |11>), which 
has Schmidt rank 4, but is less robust than the GHZ state 
of 4 qubits. 

There is one more class of states for which we were 
able to study the entanglement robustness for a general 
number of qubits n, in particular n tending towards in- 
finity, namely states exhibiting spin squeezing JToj ] . An 
n-qubit state \tp) is called spin squeezed, if it satisfies 



e = 



w a +w 9_<1 (8) 

for three appropriately chosen orthogonal directions x, 

n . 

y, and z. Here J x ,y.z = j J] a x,y,z are the total angular 

momentum operators, and (Jj) = (ip\J%\ip) etc. It was 
shown in [|| that spin squeezing is a sufficient condition 
for entanglement. Such states can e.g. be generated from 
an initial product state of the n qubits through time evo- 
lution with a Hamiltonian of the general form H = x^x- 
For sufficiently large n, there are states for which £ 2 
is very small, (J z ) = Cf with ( close to 1, and (J x ) — 
(J y ) = |l(J. The local depolarizing channels (]lf) can be 
quite easily applied in the present situation. If the state 
is written in the basis of Pauli matrices, then simply all 
Oi are multiplied by a factor s = 1 — d. One finds 

(Jz) - s(J z ) 

(4 2 ) = 5 + jE(^' ) )-J + jE^' ) ) = 

= (l- s ^ + s i(Ji), (9) 



where again s = 1 — d, and thus for the squeezing param- 
eter £3 after depolarization: 
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where £q is the original squeezing parameter. 

This means that under local depolarization the states 
definitely remain entangled for all scaling factors larger 
than 

1 , , 

Srrit = , (11) 

v/i + C 2 (i-^) 

This is equal to ^= = 0.707 in the limit of £ = 1 and 

£o = 0, which is approached for n — ► oo. This means 
that in this limit the spin-squeezing entanglement can 
stand more than 29 % of local depolarization. Note that 
(||) is a sufficient, but not a necessary, condition for en- 
tanglement, so the states may well be entangled for even 
larger values of depolarization. 

In the present work we have used local partially de- 
polarizing channels as our model for decoherence, corre- 
sponding to local and homogeneous noise. Our results 
emphasize that the robustness of multi-party entangled 
states will in general depend on the nature of the deco- 
herence process. In particular, GHZ states have turned 
out to be surprisingly robust under local depolarization. 
On the other hand, GHZ entanglement is of course frag- 
ile under some other operations, such as the loss of a 
single qubit, or the measurement of a single qubit in 
the {|0), |1)} basis. For the latter scenario it is however 
worth noting that the measurement has to be performed 
with certainty in order to destroy the entanglement. To 
show this consider the state generated from the n-qubit 
GHZ state \G n ) by performing local measurements in the 
{|0), |1)} basis on every qubit with a probability p. Cer- 
tain measurement corresponds to p = 1. This state is 
given by 

1 — (1 — nY l 

— L_^L(p®» + p»") + (i- p y»|G n )<G n |. (12) 

It is entangled for all p < 1, as can be seen from its partial 
transpose. Of course, the coefficient of the inseparable 
term decreases exponentially. 

One should certainly be careful in drawing general con- 
clusions from a single simple model. Nevertheless, our 
results show that multi-party entanglement can be sur- 
prisingly robust under decoherence, and that the robust- 
ness can even increase with the number of parties. This 
suggests that large-scale entanglement could be more 
frequent in natural systems than one might have ex- 
pected. Persistent entanglement would definitely require 
the presence of an entanglement-generating process - in- 
teractions between the subsystems, possibly in combina- 
tion with external excitation - that competes with the 
decoherence. We believe that the competition between 
entanglement generation and decoherence in multi-party 
systems is a promising and fascinating topic for future 
research. 
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